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Abstract 

We present the result of projected cross correlation analysis of AGNs and galaxies at redshifts from 
0.3 to 3.0. The Japanese Virtual Observatory (JVO) was used to obtain the Subaru Suprime-Cam images 
and UKIDSS catalog data around AGNs. We investigated 1,809 AGNs, which is about ten times larger 
a sample than that used in previous studies at redshifts larger than 0.6. We found significant excess of 
galaxies around the AGNs at redshifts from 0.3 to 1.8. For the low redshift samples {z < 0.9), we obtained 
correlation length of ro =5-6 h~^Mpc (7 = 1.8), which indicates that the AGNs at this redshift range reside 
in a similar environment around typical local galaxies. Wc also found that AGNs at higher redshift ranges 
reside in a denser environment than lower redshift AGNs; For z = 1.3 ^ 1.8 AGNs, the cross correlation 
length was measured as 11^3 h~^Mpc (7 = 1.8). Considering that our galaxies sample is based on optical 
observations with Suprime-Cam at the redshift range, it is expected that blue star-forming galaxies are 
dominant constituents of the clustering. It is successfully demonstrated that the use of the archive through 
the Virtual Observatory system can provide a powerful tool for investigating the small scale environment 
of the intermediate redshift AGNs. 

Key words: astronomical data bases: miscellaneous — galaxies: clusters: general — galaxies: nuclei 
— cosmology: large-scale structure of universe 



1. Introduction 

It has been believed that the origin of AGN activity 
is accretion of matter into a massive black hole at the 
center of a galaxy, and that almost all massive galaxies 
have such a black hole (e.g. Salpeter 1964; Lynden-Bell 
1969; Rees 1984; Kormendy & Richstone 1995; Richstone 
et al. 1998). In order to explain the activity of the AGNs 
with mass ^ 10^ M©, a large fraction of matter in the 
galaxy must be delivered to the inner region of the galaxy 
on a short timescale, ^ 10^ yr (Hopkins et al. 2008). One 
possible mechanism for causing rapid gas inflows toward 
the central region would be a major galaxy merger be- 
tween gas-rich galaxies (e.g. Sanders et al. 1988; Barnes 
& Hernquist 1991; Barnes & Hernquist 1996; Kauffmann 
and Haehnclt 2000; Granato et al. 2004; Springel et al. 
2005; Croton et al. 2006). If this is the case, AGNs arc 
expected to be found in an environment ^ of higher den- 
sity than typical galaxies since the merger occurs more 
frequently at a higher density region. 

There is some observational evidence that mass assem- 
bly of a larger system terminates earlier than that of a 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan (NAOJ). 
Data is retrieved from the JVO (http://jvo.nao.ac.jp/portal) 
operated by NAOJ. 

The term of environment is used to represent the excess number 
density around the region relative to the average at the redshift. 



smaller system (e.g. Cowie et al. 1996; Kodama et al. 
2004). From these observational facts, it is expected that 
at higher redshift the merger occurs in an environment 
with high galaxy density. Thus the AGNs produced at 
the higher redshift should be observed in a more crowded 
region. 

Measurements of the two-point correlation function of 
optically selected QSOs were performed by using the 
large dataset of QSOs obtained by 2dF and SDSS sur- 
veys (Croom et al. 2005; Ross et al. 2009; Shen et al. 
2009; Myers et al. 2006; Myers ct al. 2007). Ross et al. 
(2009) measured the real-space correlation lengths for 
z =0.3-2.2 QSO, which was in the range of 6-11 h~^Mpc 
for fixed power index of 7 = 2.0. Shen et al. (2009) an- 
alyzed data of QSOs at higher redshifts and obtained 
the real-space correlation lengths of 14.58±2.7 h~^Mpc 
(7 = 2.0) and 21.04±3.39 h-^Mpc (7 = 2.0) at z =2.9-3.5 
and z =3.5-5.0, respectively. 

The correlation length of the local quiescent galaxy 
has been measured by many authors (e.g. Brown et al. 
2000; Norberg et al. 2002; Budavari et al. 2003; Madgwick 
et al. 2003; Hawkins et al. 2003; Zehavi et al. 2005; Ma 
et al. 2009). Zehavi et al. (2005) measured real-space cor- 
relation length and obtained 5.59±0.11 h~^Mpc for flux- 
limited sample with 14.5 <r< Yl .11 and -22 < < -19 
(for h = 1). Thus the large scale environment of QSOs is 
almost equivalent to or only slightly larger than the envi- 
ronment of typical local galaxies at redshift less than 2.2, 
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and at higher redshift the galaxy number density around 
the QSOs increases with the redshift. 

The direct measurement of QSO two-point correlation 
function requires a large number of QSO samples and large 
sky coverage, thus only the large survey projects such as 
2dF and SDSS can provide a precise and reliable mea- 
surement. The measurement of QSO-galaxy cross correla- 
tion function, on the other hand, provides a direct insight 
on the local environment of QSO with fewer samples, al- 
though it is limited to lower redshifts where spectroscopic 
observation of galaxies is feasible. 

A lot of studies have been made to investigate the re- 
lation of QSO/AGN environment with a type and lumi- 
nosity of the AGNs and their evolution with redshift (e.g. 
Miller et al. 2003; Sorrentino et al. 2006; Coldwell and 
Lambas 2006; Serber et al. 2006; Strand et al. 2008; Barr 
et al. 2003; Coil et al. 2007; Adelbcrger and Steidcl 
2005; Bornancini et al. 2007). 

Studies of the environment around optically selected 
AGNs at redshift < 0.6 have been carried out by many 
authors by using the SDSS data. Miller et al. (2003) es- 
timated the fraction of galaxies associated with an AGN, 
based on 4921 galaxies with redshift less than ^ 0.1. They 
found that the fraction of galaxies associated with an 
AGN is independent of density of local galaxy. Sorrentino 
et al. (2006) examined 1829 AGNs at redshift < 0.1, and 
found no evidence of a relation between the density of lo- 
cal galaxy and the AGN activity. Coldwell and Lambas 
(2006) analyzed the environment of -^2000 AGNs at red- 
shift < 0.2, and found that the number density of galaxies 
around the AGNs is similar to that of a typical galaxy. 

Although these results show that the optically selected 
AGNs do not reside in an area of high overdensity, there 
are several studies which indicate the small-scale over- 
density around the AGNs. 

Serber et al. (2006) analyzed small-scale environment 
of 2,028 z < 0.4 SDSS QSOs and found that the QSOs 
are located in higher density region than are L* galax- 
ies and the over-density increases with decreasing scale at 
a distance of less than 0.5h~^Mpc of the QSO. They also 
reported that there is a luminosity dependence of the den- 
sity enhancement at the small-scale and the brightest end. 
Strand et al. (2008) also investigated the relationship of 
the AGN environment to type (type I, type II), luminos- 
ity and redshift of the AGN at small distance less than 
2 h~^Mpc. They extended the redshift range up-to 0.6 
compared to the work of Serber et al. (2006) by applying 
a photometric redshift cut in their selection of galaxies to 
reduce the projection effect. The number of AGN samples 
was increased to ~ 11000. They concluded that the AGN 
environment depends on the luminosity of AGNs but not 
on the type of AGNs, and that there is marginal evidence 
on variation of its environment around type I QSOs as a 
function of redshift. 

Since the observation of galaxies by the SDSS is limited 
to the redshift of less than ^^0.6, measurements of AGN 
environment at higher redshifts need to be carried out 
by an infrared instrument on a 2-4 m class telescope or 
optical and/or near infrared observations by 8-10 m class 



telescopes. 

Coil et al. (2007) measured clustering properties around 
52 SDSS and DEEP2 QSOs at redshifts 0.7 - 1.4, by 
cross-correlating each QSO with DEEP2 galaxies. They 
found that the clustering scale length between the QSOs 
and galaxies is comparable to that obtained by auto- 
correlation of DEEP2 galaxies, and obtained the corre- 
lation length of 3.3±0.7 h~^Mpc. They found no signif- 
icant dependence on luminosity nor redshift. Adelbcrger 
and Steidel (2005) examined clustering property of 79 
AGNs with absolute magnitude from —30 to —20 at 
redshifts from 1.5 to 3.5, by cross correlating with the 
1627 galaxies with measured redshift. The galaxy sam- 
ple used in their analysis was constituted of color-selected 
spectroscopic targets. They divided 79 AGNs into two 
luminosity groups, and derived correlation lengths of 
4.7±2.3 h-^Mpc and 5.4±1.2 h'^Mpc for the low and high 
luminosity groups, respectively. They reported no evi- 
dence of luminosity dependence. Bornancini et al. (2007) 
examined the correlation between the SDSS QSOs and 
distant red galaxies (DRG) at redshifts from 1 to 2, and 
obtained correlation length of 5.4±1.6 h^^Mpc. 

All of these results for high redshift (z =0.7-3.5) 
AGN/QSOs also show no significant over density of 
their environment. Although the optically selected 
AGNs/QSOs seem to be located in an environment com- 
parable to that of local typical galaxies as introduced 
above, there have been some studies which report some 
types of QSOs such as radio-loud QSOs and X-ray se- 
lected AGNs are associated with a higher density region. 

Barr et al. (2003) observed environments of 21 radio- 
loud QSOs at 0.6 < z < 1.1 with multiple filters of optical 
to near infrared bands. The observation showed that the 
radio-loud QSOs at this redshift range exist in a wide va- 
riety of environment (field galaxies, compact groups, and 
rich clusters), and indicated no evidence of redshift de- 
pendence of the environment. They also found that the 
QSOs are not always located at the center of a galaxy 
group or cluster, which contrast with low redshift clusters 
where QSOs are located at the center of galaxy distribu- 
tion. Enoki et al. (2003) predicted that, at redshift from 1 
to 2, QSOs are located in various environments from small 
groups of galaxies to clusters of galaxies based on a semi- 
analytic model of galaxy and QSO formation. This predic- 
tion is consistent with the observational evidence obtained 
by Barr et al. (2003). The over density around the radio- 
loud QSOs were also reported by some other authors (e.g. 
Yee & Green 1984; Yee & Green 1987; EUingson et al. 
1991; Hutchings et al. 1999; Kauffman et al. 2008). 

In contrast, several studies reported that the envi- 
ronment of the radio-loud and radio-quiet QSOs were 
indistinguishable (e.g. Fisher et al. 1996; Wold et al. 
2001; McLure & Dunlop 2001), though the number of 
samples used by these authors was not so large as several 
tens of samples. 

There have been a lot of studies on the environment of 
the X-ray selected AGNs (e.g. Mulhs et al. 2004; Grazian 
et al. 2004; GiUi et al. 2005; Yang et al. 2006; Miyaji et al. 
2007; Hickox et al. 2009; Gilh et al. 2009; Krumpe et al. 
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2010; Coil et al. 2009). 

Hickox et al. (2009) examined the environment of AGN 
selected from radio observations with Weterbork Synthesis 
Radio Telescope, X-ray observations of Chandra XBootes 
Survey, and mid-IR observations of Spitzer IRAC Shallow 
Survey. They measured cross-correlation between these 
AGNs and galaxies on scales of 0.3-10 h~^Mpc and de- 
rived the correlation length of 6.3±0.6 (7 = 1.8 ±0.2), 
4.7±0.3 (7 = 1.6±0.1), and 3.7±0.4 h-^Mpc (7 1.5±0.1) 
for AGNs selected in radio, X-ray, and infrared band. 
They concluded that the X-ray selected AGNs were clus- 
tered similar to typical galaxies, less clustered than radio 
AGNs, and more clustered than IR AGNs. 

Krumpe et al. (2010) examined the environment of 
X-ray selected AGNs, by using cross correlation mea- 
surement between ~ 1550 ROSAT broad-line AGN and 
~ 46,000 luminous red galaxies (LRGs) at z =0.16- 
0.36. They obtained real space correlation length of 
6.93to:28h"^Mpc for the AGNs and LRGs, and also esti- 
mated the AGN auto correlation length of 4.3]']Q gh~^Mpc. 
They also measured an X-ray luminosity dependence of 
the clustering and concluded that low luminosity sam- 
ples have a correlation length similar to that of blue star- 
forming galaxies at low redshift, and high luminosity sam- 
ples have a larger correlation length than the lower lumi- 
nosity samples and is consistent with the clustering of red 
galaxies. 

The other studies which measured auto correlation of X- 
ray selected AGNs from ROSAT all sky survey by MuUis 
et al. (2004) and Grazian et al. (2004), from Chandra Deep 
Field North and South by Gilli et al. (2005) and Yang et al. 
(2006), from XMM-COSMOS field by Gilh et al. (2009) 
and Miyaji et al. (2007), and from AEGIS field by Coil 
et al. (2009). They obtained auto-correlation lengths of 
e-lOh-^Mpc for the X-ray selected AGNs. These resuhs 
are slightly larger than typical local galaxies and quite 
similar to the optically selected SDSS QSOs. 

Although the observations of AGN-galaxies clustering 
up to redshift ~2.0 indicate that AGNs do not reside in a 
particularly high density environment, QSO-QSO cluster- 
ing study indicates that the environment of AGN at higher 
redshift is more crowded than in the low redshift Universe. 
Considering that the number density of bright AGNs 
peaks around redshift ^2 (Ueda et al. 2003; Richards et al. 
2006), it is important to extend the study of AGN environ- 
ment up to redshift 2 and beyond to solve the mechanism 
of co-evolution of AGNs and galaxies. 

In this paper we present the result on measurements of 
AGN environment at redshifts from 0.3 to 3.0, which is the 
first scientific result obtained from the Japanese Virtual 
Observatory (JVO). The reduced Subaru Suprime-Cam 
images are open to public through the JVO system. The 
JVO portal provides a functionality to search for multi- 
ple regions from its image data service by issuing a single 
JVOQL, which is an extension of database language SQL 
and can describe a coordinate join between the Suprime- 
Cam image metadata table and a QSO/AGN catalog ta- 
ble. The UKIDSS World DR2 is also used to detect galax- 



ies which arc dark in optical bands but bright enough to 
be detected in near-infrared bands. The total number of 
QSO/AGN samples used in this work is - 1809, and --90% 
of the QSO/AGN samples comes from the SDSS and 2dF 
QSO catalog. Note that this work uses the largest sam- 
ple with deep optical images, which are typically deeper 
than 24 mag. By using the deep optical images, we can 
also measure the clustering property of faint and blue star 
forming galaxies around AGNs at high redshift {z^ 1 ~ 2), 
which have not been well explored by any other studies. 
Note also that our work would be free from the cosmic 
variance, since the AGNs we have used are distributed 
through a wide of the sky. 

Throughout this paper, we assume a cosmology with 
r2,„ = 0.3, = 0.7 and h ~ 0.7. All magnitudes are given 
in the AB system. All the distances are measured in co- 
moving coordinates. The correlation length is presented 
in unit of h~^Mpc. 

2. The Datasets 

All the data used in this work were obtained 
from the Japanese Virtual Observatory (JVO) portal 
(http://jvo.nao.ac.jp/portal). The datasct downloaded 
through the JVO portal is: Catalog of Quasars and Active 
Galactic Nuclei (12th Ed.) by Veron-Cctty et al. (2006), 
SDSS DR-5 Quasar Catalog (4th Ed.) by Schneider 
et al. (2007), Subaru Suprime-Cam reduced images (ver- 
sion 0.1.24) of the JVO Subaru archive, and UKIDSS DR2 
catalog by Warren et al. (2007). Now we describe the pro- 
cedure for selecting data in order to achieve our research 
objectives. 

First we searched for AGNs observed with the Suprime- 
Cam or the UKIDSS survey at the redshift range of 0.3- 
3.0. We selected AGNs whose K-corrected absolute V 
band magnitude was in the range from —30 to —20. 
The K correction was made by assuming power law in- 
dex of a = —0.5 for AGN spectrum and calculated as 
k = -2.5(a -I- 1) log(l -I- z). Then we removed AGNs that 
were located near another lower redshift AGN within an 
angular distance corresponding to 8 Mpc at the AGNs' 
redshift. This filtering aimed to reduce the projection 
effect of galaxy clustering associated with other nearby 
AGNs. 

For each AGN among these samples, the Suprime-Cam 
images of B, V, Rc, Ic, i', and z' band and the UKIDSS 
K-band data were retrieved, and uniformity of the images 
and data coverage for the AGN field were examined. 

The quality of some of the Suprime-Cam images re- 
trieved from the JVO is highly inhomogeneous, because 
the images were created by stacking all the frames. 
Certain parts of an image, therefore, may have signifi- 
cantly shorter exposure time than other parts. The in- 
homogeneity was estimated from the distribution of root 
mean squares (RMS) of the background noise, which was 
calculated by the SExtractor (Berlin and Arnouts 1996) 
at every image pixel. We removed images where the stan- 
dard deviation of RMS over the entire AGN field was 
larger than 10% of the average. The Suprime-Cam im- 
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ages which were not associated with photometric cahbra- 
tion data were not used in this analysis either. 

Since we were interested especially in an environment 
of high redshift (z > 1) AGN, we masked bright and ex- 
tended sources on the Suprime-Cam images, so that low 
redshift bright galaxies were effectively filtered from our 
galaxy sample. The masked regions were defined as a re- 
gion where more than 100 contiguous pixels were illumi- 
nated with a flux higher than 10 sigma of the background 
fluctuation. Comparing with an unmasked image, we es- 
timated that the masked sources were about two magni- 
tudes brighter than the peak of the magnitude distribu- 
tion, which was ~22 mag on average. Object catalogs for 
Suprime-Cam images were created separately for each ob- 
servation band using the SExtractor (Bertin and Arnouts 
1996). 

In order to examine the data coverage of the AGN field, 
an effective observation area was calculated for each annu- 
lus centered on an AGN location, which has a correspond- 
ing width of 0.1 Mpc. The effective area observed with the 
Suprime-Cam was calculated by counting the unmasked 
pixels of the image. The coverage of the UKIDSS data 
was estimated from a fraction of an area where number 
density of objects was lower than the average at the same 
annulus by five sigma. 

It is difficult, however, to measure the coverage fraction 
of the annulus for which the expected source number is 
too small to measure the density variation in the annulus. 
To evaluate the coverage of the low count annulus, we es- 
timated the average number density expected for full cov- 
erage, and compared the observed number in the annulus 
with the expected number. The average number density 
was estimated from the annulus with second largest den- 
sity among the top-five annuli in terms of the number of 
detected sources. If the observed number is below the ex- 
pected number by five sigma, we considered the annulus 
to be an unobserved area, that is 0% of coverage. 

We set a threshold on the data where at least 80% of 
every annulus should be observed. The data that did not 
meet this threshold were removed. The fraction of data 
coverage estimated as described above was used in the 
calculation of number density of galaxy at each annulus. 
After these selections, 2689 AGNs were selected for further 
analysis. 

To improve the signal to noise ratio, we removed the 
data of shallow observations based on the average number 
density, poj of detectable galaxies at the AGN redshift. po 
was estimated by integrating a luminosity function up to 
an effective limiting magnitude, mumit, of each Suprime- 
Cam image or UKIDSS data. The detailed description for 
determination of po will be described in the next section. 

Then we selected data where po was greater than 
10~^ Mpc~^ at any one of the observation bands. The 
distribution of po vs AGN redshift z are shown in fig- 
ure 1 for data of po > 10~^. In the figure, po is the largest 
value among all the observation bands for each AGN. One 
can see a clear separation between the samples which are 
observed in the UKIDSS Large Area Survey (LAS) (lower 
cluster of open circles), and the samples observed with the 



Suprime-Cam and UKIDSS Deep Extragalactic Survey 
(DES) (upper cluster of crosses and open circles). We 
applied the selection criterion of po > 10~^ Mpc""^, which 
limit the maximum redshift of UKIDSS LAS samples to 
^1.3. The number of AGN selected with this criterion is 
2023. 

The foreground cluster or group of galaxies can be a 
noise for the clustering analysis. To reduce the effect of 
accidental alignment of the nearby cluster, we calculated 
a clustering coefficient, Bqq, around each AGN, which 
was defined as ^(r) = Bqgt'^ , and Bqg was calculated 
as (Barr et al. 2003): 

Bqo = i-:2 ^*°^"'-^^^ (lMpc)-^-^ (1) 
27r Po 

where Atotai is the total number of observed galaxies at 

< 1 Mpc, and Nj^g is the expected background count 
at Tp < 1 Mpc estimated from the averaged density at 
fbg,min <rp < rbg,max- The distance ranges rbg,min and 
?'bg,min wcrc determined for each redshift group, and they 
are summarized in table 1. If a cluster of galaxies is 
accidentally located in front of the AGN, significantly 
large Bqg values would be observed. Figure 2 shows the 
distribution of log | Bqg I for 2023 AGNs which satisfied 
the Po criterion. The selection criteria used here were 
\Bqg\ < 1 X 10^. Abell class clusters corresponds to 
Bqg 300 Mpc^-^^ (McLure & Dunlop 2001). The up- 
per bound of this selection range corresponds to a cluster 
which has 30 times larger clustering coefficient than the 
Abell class 0, and it is expected that Bqg of a real system 
rarely exceeds this criteria. Thus most of the samples re- 
jected with this criteria are the ones with an uncertainty of 
the clustering coefficient as large as this boundary or clus- 
ters accidentally located near the line of sight. Applying 
both the Po and Bqg criteria, the number of AGNs was 
reduced to 1976. 

It is also possible that the nearby cluster or high den- 
sity region is located in a region offset from the AGN. To 
remove such samples, the reduced and maximum devi- 
ation Gmax of radial distribution of galaxy density around 
the AGN were calculated at distance range of 1 Mpc to 
''bg.min- The distributions of the reduced and maximum 
deviation Cmax for the 1976 AGN samples are shown in fig- 
ures 3 and 4. The criteria used here were x^/(n — 1) < 3.0 
and (Jmax < 5.0. After this selection, the total number of 
AGNs was 1828. 

Figure 5 shows the distribution of K-corrccted V-band 
absolute magnitude and redshift of the AGNs selected by 
the above criteria. Open circles represent AGNs for which 
the analysis was made using the UKIDSS data, and crosses 
represent AGNs for which the analysis was made using the 
Suprime-Cam data. The histograms of AGN redshift for 
dim AGNs (My > -25) and bright AGNs {My < -25) are 
shown in figure 6. 

In table 2, the number of AGNs used in this work are 
shown for each data subset (zl-D to z5-B defined in sec- 
tion 4) and for each origin of the AGN identification. The 
low redshift and high luminosity AGNs, and high redshift 
and low luminosity AGNs were not analyzed for their poor 
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statistics. As a result the total number of analyzed AGNs 
was reduced to 1809. 

The column "SDSS" of the table represents the number 
of AGNs which are contained in the SDSS QSO catalog 
(4-th ed.). The column "2dF" represents the number of 
AGNs which are contained in the 2dF QSO catalog but 
not contained in the SDSS catalog. The AGNs which are 
contained in both the SDSS and 2dF catalog are counted 
as "SDSS" . Among the AGNs which are not contained in 
the SDSS catalog nor the 2dF catalog, the AGNs which 
were identified based on an optical, UV, or infrared ob- 
servations are counted as "UV-OPT-IR" , the AGNs which 
were selected from X-ray or Radio catalog are counted as 
"XRAY" or "RADIO", respectively. About ninety per- 
cent of the AGNs are selected from the large survey of 
SDSS and 2dF, thus our AGN samples are nearly homo- 
geneous and dominated by optically selected AGNs. 

The celestial distribution of the 1809 AGNs is shown in 
figure 7. Since a large fraction of AGNs comes from the 
SDSS catalog, the distribution almost follows the SDSS 
survey area. The distribution of the right ascension and 
the declination are also shown in figure 8 for all the sam- 
ples (open histogram) and samples which have data of 
Suprimc-Cam (shading histogram). The concentration 
to declination ^ degree is due to the contribution of 
UKIDSS LAS survey. The most significant concentration 
of Suprime-Cam observations is located at the Sextans 
field (RA= 150° - 157°, Dec= -4° - -f 1°). The total num- 
ber of AGNs in the Sextans field is 80, and their number 
fraction at each subset group is ~ 20% for z4-B and z5-B 
group, and less than 6% for the other groups. We con- 
firmed that this concentration to a specific field did not 
affect the result so much, and it will be discussed in sec- 
tion 5. Although there is a slight concentration to specific 
fields, the AGNs are distributed in an extended area of 
the sky, and we can expect that the bias effect (cosmic 
variance) would be mostly negligible in our analysis. 

The Suprime-Cam images provided by the JVO por- 
tal were processed using the JVO Suprime-Cam reduction 
pipeline (version 0.1.24) which was developed based on the 
Suprime-cam Deep Field REDuction (SDFRED) package 
(Ouchi et al. 2004). The data reduction was carried out 
on the JVO grid computing system (Shirasaki et al. 2008). 
All the data obtained before December 2006 were used to 
create mosaic images. The SDFRED is based on a soft- 
ware package developed for Suprime-Cam data reduction 
(Yagi et al. 2002). 

The data were reduced in the following procedure for 
each CCD frame: 1) bias subtraction; 2) fiat field correc- 
tion; 3) distortion correction; and 4) astrometric correc- 
tion. The positions, rotation angles, and flux normaliza- 
tion factors relative to a reference frame were calculated 
for every frame, then they were stacked to make a single 
mosaic image. The photometry was performed using the 
data of standard star observations which were taken on 
the day that one of the stacked frames were observed. The 
comparison with the SDSS catalog shows the differences 
of photometric magnitudes were less than 0.2 for 90% of 
the data compared. The astrometric accuracy was esti- 



mated by comparing with the USNO catalog. The peak 
of the distribution of angular difference from the USNO 
data was 0.3 arcsec, and the difference was less than 1.2 
arcsecond for 99% of the analysed data. 

3. Estimation of po 

Estimation of po, which is an average number density 
of observed galaxy at AGN redshift, is crucial for deter- 
mining the cross correlation length of AGNs and galaxies. 
Since redshift of galaxy is not measured in this work, po 
cannot be determined directly from the data used here. 
Therefore, it needs to be indirectly derived by model cal- 
culation under reasonable assumptions. 

The luminosity function of galaxy (t>{M) is one of the 
functions required for the model calculation. Only the ob- 
servable that can be used to infer the density of observed 
galaxy is an observed magnitude distribution N{m). The 
observed magnitude distribution is a product of true 
magnitude distribution 7Vtruc(w) and detection efficiency 
DE{m) for an object of apparent brightness m. Once the 
model functions for (l){M), A'truc("i) and DE{m) are de- 
rived, we can estimate po by comparing the observable 
N{m) with the model prediction. 

A brief procedure for determining po and its uncertainty 
is described as follows; we define a parameter mumit which 
is closely related with the limiting magnitude of observa- 
tion and is used to calculate po by integrating the lumi- 
nosity function (j){M) to the corresponding absolute mag- 
nitude at the AGN redshift. We estimate mnmit under 
the assumption than the cumulative fraction of apparent 
magnitude above mpcak has a constant value of i^umit at 
wiimit- This assumption would be valid if the observations 
were carried out under identical conditions. This is, how- 
ever, not actually the case, and Fumit is not a constant but 
a variable depending on the observational condition. Thus 
we calculated the expected value of Fiimit for each obser- 
vation of AGN field, and we considered the most frequent 
value of -Fiimit to be an estimator for the most probable 
value of Po, and the lower and upper limit to i^umit to be 
an estimator for uncertainty of po. 

First we describe the model of luminosity function used 
in this work. We parametrized the luminosity function of 
galaxy by means of the absolute magnitude M, redshift 
z, and source frame wavelength A by using the luminosity 
functions obtained by Gabasch et al (2004) for 1500A, 
2800A, u, B and g' bands, Gabasch et al (2006) for r', i', 
and z' bands, and Cirasuolo et al. (2007) for K band. The 
parameters a, (j)o, and Mq of the Schechter function for 
each rest-frame wavelength band were parametrized by a 
polynomial function of redshift. Then, for each redshift, 
these parameters were represented as a function of rest- 
frame wavelength by interpolating with the cubic spline. 
Using this parametrization, we obtained luminosity func- 
tion 4>{M] z, A) for an arbitrary redshift and rest-frame 
wavelength. 

Po is calculated by integrating the product of the 
luminosity function 4>{M; z, A) and detection efficiency 
DE{M + DAI) for a source with absolute brightness of 
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M (apparent brightness m = M + DM) at AGN redshift 
z. DM is a distance modulus for redshift z. For sim- 
plicity, we use an effective limiting magnitude miimit as a 
characteristic parameter for DE{m). mii,„it is defined so 
that the integral of luminosity function up to mumit ~ DM 
equals to pq: 

^-DM 



Po = 



(j3{M)DE{M + DM)dM^^ 



low ■ 



-DM 



^limit 



'low 



-DM 



where miow is a lower boundary of the apparent magnitude 
for the integration. 

The luminosity functions were evaluated for the AGN 
redshift and for the rest frame wavelength corresponding 
to the center wavelength of the observation band. We 
assumed the following function for modeling the detection 
efficiency : 



DE{m) = 



1 



Cxp(-(TO-mth)^/cr™) 



(to < mth) 
(to > mth), 



(3) 



where mth and dm represent the magnitude where de- 
tection efficiency starts to decrease and the attenuation 
width of the detection efficiency, respectively. In the case 
of Suprime-Cam data, the detection efficiency rapidly de- 
creases at the bright magnitude side owing to the masking 
procedure described in the previous section. This effect 
was approximated by introducing the lower cut-off magni- 
tude TOiow- ni\ovi is a minimum of the observed magnitude 
excluding the first 5% of all the sources ordered with their 
magnitudes. The number density at the bright end of the 
luminosity function is small and the error on the po es- 
timation is expected to be lower than a few %, which is 
negligible compared to the error caused from the uncer- 
tainty of mnmit- 

The only observable that can be used to determine 
wiiinit is the observed magnitude distribution N{m). 
Especially, the distribution beyond the peak magnitude 
Wpeak can be a sensitive estimator for mhmit- Thus we 
define the cumulative fraction F{m) to estimate TOumit as 
follows: 



F{m) = / 

J rn 



N{m')dm'/ 



peak 



N{m')dm'. 



(4) 



are shown in the figure. The best fit parameters for 
these examples are: a = 0.47 ± 0.02, b = 0.20 ± O.OI, 
TOb = 19.0 ± 0.1, mth = 20.8 ± 0.06 and cr„ = 1.11 ± 0.03 
for UKIDSS K-band, and a = 0.98 ± 0.05, b = 0.28 ± 0.01, 
mb = 23.4 ± 0.04, mth = 26.1 ± 0.02 and ct™ = 0.51 ± 0.01 
for Suprime-Cam R-band. The break magnitude mb ap- 
proximately corresponds to Mq of the Schechter function 
(j}{A§)dM'^^^™^'i^ observable redshift, thus it is expected 
to correlate with mth which is a parameter related with 
sensitivity of the observation. Since an accurate value of 
mb is not required for the estimation of mhmit, we simply 
used an empirical relation of mb — mth — 2.5. 

The cumulative fraction at TOumit, ^iimit = ^'(mhmit), 
depends on the limiting magnitude of the observation, red- 
shift of the AGN, the observation band, and the model 
parameter of A^truo- We, therefore, calculated Fumit for 
various plausible parameter ranges: z = 0.3-3.0, a„i = 0.5- 
1.2, a ^ 0.6, 6 = 0.1 - 0.6, mth = 19 - 26, and po > 10"^. 
The last criteria is the same as was applied to the obser- 
vational data. To choose the ranges of parameters (Tm, b 
and mth, we fit the magnitude distribution for each AGN 
sample to the model function iVmodcKm). Since the statis- 
tics of UKIDSS LAS samples were too poor to fit to the 
model function, we summed every five samples in order 
from the lowest mpcak sample, and fit the model to the 
summed distribution. The ranges of 6, am and mth were 
determined so that more than 90% of the samples were 
included. Since the parameter a is insensitive to i^hmit, it 
was fixed to an average value of 0.6. 

For these parameter ranges, we obtained 0.1, 0.6, and 
0.8 for the minimum, median, and maximum values of 
-Fiimit, respectively. Therefore, we calculated the most 
probable value of po by integrating the luminosity function 
up to TOiiniit determined with i^iimit — 0.6, and lower and 
upper boundary with -Fiimit = 0.1 and 0.8, respectively. 

4. Cross-Correlation Analysis 

The two-point cross-correlation function, ^, is gener- 
ally used to measure the clustering properties between an 
AGN and its surrounding galaxies. It is defined as: 



'peak 



To calculate a model function for F{m), we need a model 
function for true magnitude distribution A'truc(m). We 
assumed a broken power law function for Ntmc{m): 

c-10''("-™b) (m<mb) 
c-10''(™-™b) (m>mb). 



(5) 



where TOb is a break magnitude of the magnitude distri- 
bution. Using this model function, the observed mag- 
nitude distribution can be calculated as A^modci(m) = 
A^truG(m) X DE{m). Applying this to the equation (4), 
we can obtain the model function -Fmodei(™)- 

To show how well the model function A^modci(m) de- 
scribes an actual observed magnitude distribution, the 
model functions fitted to the observation are shown in 
figure 9. The fitting results for two examples of UKIDSS 
K-band (left histogram) and Suprime-Cam R-band (right) 



^(r)=p(r)/po-l, 



(6) 



where p{r) is the number density of galaxies at a distance 
r from an AGN, and po is the average number density 
of galaxies at the AGN redshift. Since the redshifts were 
not measured for galaxies, the projected cross-correlation 
function, w, was derived from the observational data in 
the following way. 

The projected cross-correlation function is obtained by 
integrating the equation (6) along the line of sight, then 
it is described as (Davis & Peebles 1983): 

/>oo />oo 

u{rp)^2 C(rp,7r)d^ = 2 / rdr^(r)(r2 -r2)-"-5,(7) 

Jo Jvp 

where rp and tt is components of distance from the AGN 
perpendicular to the line of sight and along the line of 
sight, respectively. We assumed a power law function for 
the density profile around an AGN, e.g. ^(r) = (ro/r)'^. 
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In this case, equation (7) can be integrated analytically 
as: 

V„\^r(i/2)r((7~i)/2) 



(8) 



'P/ r(7/2) 

where r is the Gamma function. 

Studies of cross-correlation analysis published usually 
used a method which counted pairs of objects with a 
given separation {DD{r)) and compared with the ex- 
pected background which is a number of pairs in a random 
sample for the same separation {DR{r)) (Davis & Peebles 
1983), 

DD{r) 



C(r) = 



1. 



(9) 



DR{r) 

In this work, however, data of each AGN field are highly 
inhomogeneous, thus it is not appropriate to use the ran- 
dom catalog created from the whole dataset for estimating 
the background level. In addition, we need to consider two 
kinds of background. One is a background corresponding 
to DR{r) in the equation (9), which is an average num- 
ber density at the AGN redshift. We estimate this back- 
ground po as explained in section 3. The other background 
is an integrated number density of all the observed galax- 
ies from z = to infinity, which needs to be considered 
in our case as the redshifts of galaxies were not measured 
in our sample. The second background, ribg, is estimated 
by fitting model function to the observed data. In ad- 
dition, the coordinates of AGNs in our samples are well 
separated from each other, and the correlation function 
for each AGN can be calculated independently from the 
other AGN samples. Thus, in our case, it is enough to 
simply calculate the number density around each AGNs 
and take an average of them. 

We measured number density of all the detected sources 
around AGNs of a given redshift range. The number den- 
sity at perpendicular distance rp from the i-th AGN is 
calculated as: 

n,irp)=N,{rp)/S,{rp), (10) 

where Ni{rp) is the number of detected sources at a dis- 
tance range of rp^ rp^ dr for the i-th AGN, Si{rp) is the 
projected area corresponding to the distance range. Then 
the number density averaged over m AGNs of the same 
redshift range are given by: 



n(rp) 



— Vn,(rp) 

771 ^ ^ 



(11) 



The projected correlation function w(rp) averaged over 
the m AGNs was calculated as: 



uj{rp) 



n{rp) - Hbg 



(12) 



Pi^o is the averaged number density of observed galax- 
ies at the i-th AGN redshift. pi^ was calculated for the 
rest frame wavelength corresponding to each observation 
band as described in section 3. The rest frame wavelength 
corresponding to the observation band was calculated by 
Arcst = Acontor/ (1 + z), whcre Acontcr IS thc effective center 
wavelength of the observation band. 



We used only the data observed at the most sensitive 
band for detecting the galaxies at the AGN redshift. The 
observation band with the largest value of p^.o was chosen 
as thc most sensitive band. To show the depth of obser- 
vations used in this work, miimit of the observation with 
the most sensitive band is plotted against redshift in fig- 
ure 10. Open circles represent AGNs for which UKIDSS 
survey data were used, and crosses represent AGNs for 
which Suprime-Cam data were used. The apparent mag- 
nitude corresponding to an absolute magnitude of —21 is 
shown as a dashed line. The number of AGN samples for 
each observation band arc summarized in table 3. 

In calculating Ni{rp), we took into account the coverage 
fraction for each annulus, fi{rp), by multiplying 1/ fiirp) 
to a number of objects detected at a distance bin Vp ^ 
Tp + dr. 

Galaxies near an AGN may be contaminated by the 
light of the AGN, and the detection efficiency for galax- 
ies might be lowered there. Thus, in calculating Ni{rp), 
the region of 4 arcsec from the AGN was masked. The 
typical FWHM of thc Suprime-Cam image used in this 
work is less than 1 arcsec, therefore the contamination ef- 
fect of the AGN light to the galaxy detection efficiency 
would be negligible outside the masked region if the ap- 
parent brightness of the AGN is dimmer than 17 mag. 
The reduction of thc projected area due to the masking 
was taken into account in thc calculation of Si{rp). 

5. Results 

We divided all the AGN samples into five redshift 
groups, z = 0.3-0.6 (referred to as zl), 0.6-0.9 (z2), 0.9- 
1.3 (z3), 1.3-1.8 (z4) and 1.8-3.0 (z5). For each red- 
shift group, the dataset is further divided into a dimmer 
group {My > —25, referred by putting a suffix D after 
the redshift group name, e.g. zl-D) and a brighter group 
(My < —25, suffix B). In figure 5, the redshift and luminos- 
ity range of each group is shown with a solid box together 
with the absolute magnitude and redshift distribution of 
the AGNs. Hereafter we call an AGN of My > -25 as 
"dim AGN", and an AGN of My < -25 "bright AGN". 
The bright AGNs for z < 0.6 and the dim AGNs for z 
> 1.3 were not analyzed due to their poor statistics. 

For each group the galaxy number density, n{rp), was 
calculated as a function of perpendicular comoving dis- 
tance from the AGN, and it was fitted with a model func- 
tion as described in thc previous section. In the left panels 
of figures 11 and 12 the galaxy number density for each 
AGN group is shown with closed circles, and the model 
function fitted to the observation is shown with a solid 
line. The error bars represent one sigma statistical errors. 
The model function is derived from equations (8) and (12) 
as follows: 



n{rp) : 



ro\^r(l/2)r((7-l)/2) 



Pa + ribg. 



(13) 



r(7/2) 

The horizontal dashed line in figures 11 and 12 represents 
the fitting parameter ribg of equation (13). The bin size 
was fixed to 0.1 Mpc for redshift groups of zl and z2, and 
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it was fixed to 0.2 Mpc starting from 0.1 Mpc for z3, z4 
and z5 groups. For z3, z4 and z5 groups, the data of < 
0.1 Mpc was ignored, since a large fraction of this distance 
range was masked to avoid the effect of AGN light. 

In the right panels of the figures, the corresponding 
projected correlation functions w(rp) are also shown. We 
fixed the slope parameter 7 to 1.8, which is a canonical 
value obtained from auto-correlation analysis among typ- 
ical local galaxies by many other works (e.g. Zehavi et al. 
2005). po is an average of pi^ for each AGN in the group, 
and it is fixed to the value calculated as described in sec- 
tion 3. The fitting parameters are summarized in table 4. 

Evidence of galaxy clustering around the AGNs was 
detected with a 90% confidence level (C.L.) for all the 
groups except for z5-B. We obtained correlation lengths 
of 4.?;^:^, b.StIi and T.GlfJ h'^ Mpc for dim AGN 
groups of zl-D, z2-D, and z3-D, respectively. For bright 
AGN groups, correlation lengths of 6.3l^;5, 
11. 1^2:7 h"^ Mpc were obtained for z2-B, z3-B, and z4-B, 
respectively. The errors quoted here indicate the sum of 
a systematic error estimated from uncertainty of po and 
a statistical error of one sigma. For z5-B AGN group, 
no significant clustering signature was observed, and an 
upper limit of 13.1 h^^ Mpc was obtained. This upper 
limits includes systematic error due to uncertainty of po 
and statistical error of one sigma. 

To examine the possibility that the clustering observed 
for z4-B was due to accidental alignment of high density 
regions, we analyzed eight independent offset fields of z4-B 
AGNs. The offset fields were taken apart from each AGN 
by 10 Mpc in the directions of right ascension and/or dec- 
lination, so that the offset fields did not overlap with the 
high density region around AGN. We applied the same 
analysis procedure to these offset fields as was applied to 
the AGN fields. We derived density distributions aver- 
aged over each of the eight offset fields. The significance 
of excess or deficit at r =0-2.5 Mpc to the average den- 
sity at r =4-6 Mpc were less than two sigmas, while the 
significance for the z4-B AGNs was 5.8 a. As no signif- 
icant systematic excess was seen for the offset fields, we 
conclude that the excess seen for the AGN dataset is most 
likely due to galaxies associated with the AGNs. 

We also examined the possibility that the excess was 
caused by only a single or a few samples which had large 
clustering. We calculated the average number density 
around z4-B AGNs for which absolute value of Bqq was 
less than or equal to 1000. The -Bqg distribution for z4-B 
AGN samples is shown in figure 13, and 110 out of the 
142 AGN samples satisfied this condition. The number 
density distribution is shown in figure 14. The excess at 
r < 2.5 Mpc is clearly found also for this relatively small 
clustering samples of |-Bqg| < 1000. The calculated cor- 
relation length was 8.7^2 ^^"^Mpc. From this result, we 
conclude that the excess seen in the z4-B sample is not cre- 
ated just from a single or a few samples with the largest 
clustering, but from a wide range of samples of z4-B group. 

As described in section 2, 20% of the z4-B AGN samples 
are located in the Sextans field of RA = 150° - 157° and 
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Dec = —4° - 4-1°. To see the effect of the concentration to 
this specific field, we derived the distribution of number 
density for z4-B AGNs excluding the samples located in 
the Sextans field, and showed it in the left panel of fig- 
ure 15. In the right panel of the same figure, the density 
distribution for z4-B AGNs located in the Sextans field is 
also shown. No significant difference is seen for the two 
density profiles. The cross correlation length for the z4-B 
samples excluding the AGNs of Sextans filed is lO.Olj'g 
h~^Mpc, and it is consistent with the result for the whole 
z4-B samples. Thus, we can conclude that the result of 
z4-B is not strongly biased by the samples of the Sextans 
field. 

In figure 16, the correlation length obtained for each 
AGN group is plotted as a function of redshift. The bright 
AGN samples (z2-B, z3-B, z4-B) are shown as closed cir- 
cles, and the dim AGN samples (zl-D, z2-D, z3-D) are 
shown as closed triangles. The upper limit is shown for 
the z5-B sample. For comparisons, results of AGN-galaxy 
correlation analysis by other authors (Bornancini et al. 
2007; Coil et al. 2007; Norman et al. 2009; Hickox et al. 
2009) are also shown in the figure. The auto-correlation 
length of galaxy or QSO obtained by other authors (Ma 
et al. 2009; Zehavi et al. 2005; Hawkins et al. 2003; Wake 
et al. 2008; Ross et al. 2009; Shen et al. 2007) are also 
shown in the same figure. 

It may not be appropriate to compare the clustering 
strength only with the the correlation length, since it de- 
pends on the assumed or fitted slope parameter 7. The 
clustering strength is usually compared in erg , which is the 
rms fiuctuation of the galaxy distribution in the sphere 
with a comoving radius of rmax = 8 h^^Mpc, and calcu- 
lated as (Miyaji et al. 2007; Peebles 1980). 

(t'I - (ro/w)V2 (14) 

J2 = 72/[(3-7)(4-7)(6-7)2l. (15) 

Using these formulae, we converted a correlation length of 
each experiment to ag, and they are compared in figure 17. 

The results for low redshift samples (z =0.3-0.9, zl 
and z2) are consistent with the existing measurements. 
Norman et al. (2009) measured the projected cross corre- 
lation between 420 QSOs and 4975 luminous red galaxies 
(LRG) at redshifts from 0.2 to 0.8 (an open square in fig- 
ures 16 and 17) based on the data of 2QZ and 2SLAQ sur- 
vey. Our measurements of a% &t z = 0.3-0.9 (zl-D, z2-D, 
z3-B) agree with their value, although our galaxy sample 
is a mixture of dim galaxies detected with the Suprime- 
Cam and bright galaxies detected with the UKIDSS sur- 
vey. Hickox et al. (2009) also measured the cross cor- 
relation of three types of AGNs and galaxies (asterisks 
in figures 16 and 17; from top to bottom radio. X-ray 
and IR selected AGNs) based on the AGN and Galaxy 
Evolution Survey (AGES) and Bootes multi-wavelength 
survey. Their sample contains 598 AGNs at redshift of 
0.25 - 0.8. Our results are nearly consistent with their 
result. 

The clustering of local galaxies measured by Ma et al. 
(2009) (K-band selected galaxy), Zehavi et al. (2005) (r- 
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band selected galaxy), and Hawkins et al. (2003) (6j-band 
selected galaxy) are 5-7 Mpc (open diamonds in figure 16 
from top to bottom, respectively). Our values for AGN- 
galaxy correlation length at z ~ 0.3-0.9 are, therefore, 
almost the same as the auto-correlation length of the local 
galaxies. 

From the observational facts that at the local Universe 
red early-type galaxies are more clustered than blue late- 
type galaxies and bright galaxies are more clustered than 
dim galaxies, we can expect different clustering strength 
for UKIDSS and Suprime-Cam galaxies samples. To see 
the difference, we calculated the AGN-galaxy cross corre- 
lation function for UKIDSS selected galaxies and Suprimc- 
Cam selected galaxies independently for zl-D and z2-D 
AGN groups. The results are summarized in table 3 and 
labeled as "OPT" for the Suprime-Cam galaxy samples 
and "IR" for the UKIDSS galaxy samples. 

For the zl-D AGN group, larger correlation length 
was obtained for UKIDSS galaxies samples (ro = 6.8li['o 
h~^Mpc ) than that for Suprime-Cam galaxies samples 
(ro = ^-Oti.l h-iMpc). The clustering of UKIDSS galax- 
ies around zl-D AGNs is consistent with that of K-band 
selected galaxies measured by Ma et al. (2009) at the lo- 
cal Universe. As for the Suprime-Cam galaxies, the clus- 
tering is as small as that of dim galaxies (Af^ > —19) 
measured by Zehavi et al. (2005). At this redshift range 
(zl), galaxies in the Suprime-Cam sample are dimmer 
than M ~ — 20 as the brighter galaxies were masked be- 
fore the source extraction, while the UKIDSS galaxy sam- 
ple consists of bright galaxies (M < —21) (see top panel 
of figure 18). Thus the Suprime-Cam galaxies sample is 
strongly biased to dim galaxies (expected to be in the 

range of M = -16 20), while the UKIDSS sample 

is biased to bright galaxies (M < —21). The correlation 
length obtained for the Suprime-Cam galaxies can be con- 
sistent with the optically selected dim galaxies at the lo- 
cal Universe, if the luminosity dependence measured by 
Zehavi et al. (2005) is taken into account. 

As for the z2-D group, correlation length of 7.4^2'^ 
and 5.0l};gh-iMpc were obtained for the UKIDSS and 
Suprime-Cam samples, respectively. The difference was 
smaller than that of the zl-D group, and they were nearly 
consistent within their uncertainties. At this redshift 
range, maximum brightness of the Suprime-Cam galax- 
ies comes down to M ^ — 22, which almost equals or 
is brighter than the characteristic luminosity (break 
point of luminosity function). Thus, it is expected that 
the clustering of Suprime-Cam galaxies is enhanced due 
to the increase in the fraction of bright galaxies in the 
sample. We conclude that the difference between the clus- 
tering properties of UKIDSS and Suprime-Cam galaxies 
can be explained by the luminosity dependence of galaxy 
clustering as observed in the local Universe. 

The result for z3-D is significantly larger than the result 
by Coil et al. (2007). They measured quasar-galaxy cross 
correlation length from 52 QSOs at redshifts from 0.7 to 
1.4 based on the SDSS and DEEP2 survey (an open trian- 
gle in figures 16 and 17). The range of absolute magnitude 



of AGNs used in Coil et al. (2007) is from -26 to -20, 
which is almost equivalent to our dim AGN sample. Coil 
et al. (2007) measured the correlation length using the 
data at a distance scale from 0.1 to 10 Mpc. They used 
7 as a free parameter, then obtained 1.55 for 7. Thus we 
fit the model function to our data of 0.1 to 3.5 Mpc with 
a fixed 7 of 1.55. Then we obtained much longer correla- 
tion length, which is not consistent with the result by Coil 
et al. (2007). The difference may come from the difference 
of galaxy samples used. The DEEP2 galaxies sample is 
selected by adapting a color cut, and it is expected that 
redder galaxies are more dominant in their sample than 
in ours. 

It should be noted that our clustering measurement is 
dominated by dim galaxies near the detection threshold, 
while the other experiment conducting spectroscopic ob- 
servations is mostly targeted to the bright end of the lumi- 
nosity function. Therefore, our measurement has sensitiv- 
ity to the clustering properties of relatively dim galaxies 
which are not measured by spectroscopic observations. In 
addition to that, our optical measurement is biased to 
blue galaxy at this redshift range, thus it is expected that 
star-forming galaxies are dominant in our galaxy sample. 

In the case of spectroscopic observation, a red galaxy is 
preferentially selected as a target since the redshift esti- 
mate for a red galaxy is more robust than that for a blue 
galaxy owning to the existence of a sharp 4,000 A break 
in the spectrum. 

The difference between our measurement and that of 
Coil et al. (2007) may be attributed to the difference in 
the clustering properties of a dimmer-blue galaxy and a 
bright-red galaxy. It is also possible that the difference is 
due to the cosmic variance, as the number of well sepa- 
rated fields is only four for the DEEP2 dataset. 

The cTg, AGN-G obtained for the z3-D is almost compa- 
rable with (78, Q derived from the QSO auto-correlation 
length by Ross et al. (2009), and with crg^G derived from 
the auto-correlation length of radio-quiet Luminous Red 
Galaxies (LRG) at lower redshift by Wake et al. (2008). 
The correlation length for the bright AGN group at the 
same redshift (z3-B) was consistent with that for z3-D 
group, and no significant evidence of AGN luminosity de- 
pendence was found. Shen et al. (2009) examined lumi- 
nosity dependence of the clustering of SDSS QSOs. They 
also found no significant evidence for the luminosity de- 
pendence, which is consistent with our observation. 

The distribution of galaxy number density for the bright 
AGN group at redshift 1.3-1.8 (z4-B) shows relatively flat 
excess at < 2.5 Mpc. This feature may indicate that a 
significant number of the AGNs are not located at the 
center of a galaxy group but distributed over a scale of 
5 Mpc. 

Bornancini et al. (2007) have examined the correlation 
between the 13 SDSS QSOs and the distant red galaxies 
(DRG; J— Kg > 2.3) at redshifts from 1 to 2 (an open circle 
in figures 16 and 17). Their AGN sample has a luminosity 
range similar to our bright AGN sample. They reported 
a correlation length of 5.4±1.6 h^^Mpc. Considering that 
the redshift range of their sample extends over z3 and 
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z4 of our samples, we calculated a correlation length and 
CT8,AGN-G for a combined sample of z3-B and z4-B to be 
7.5ti;gh"iMpc and LS^Qj, respectively. Although this 
result is almost consistent with the value of Bornancini 
et al. (2007) within the quoted error, the distribution of 
number density shows a flat excess at 1-2.5 Mpc (fig- 
ure 19), which was significant for z4-B sample and barely 
seen also in the z3-B sample. This feature was not ob- 
served in the data of Bornancini et al. (2007). The mag- 
nitude range of galaxy sample used by Bornancini et al. 
(2007) is K < 19.5(Vega) - 21.4(AB), while our galaxy 
sample of Suprime-Cam is dimmer than mAB ~ 22 on 
average. Thus the difference in the number density pro- 
file might be due to the difference in brightness and/or 
observation band of the galaxy samples. 

The clustering strength, cts^agn-g, obtained for z4-B 
sample is nearly consistent within the quoted error with 
the QSO clustering crg^q calculated from the observational 
result by Ross et al. 2009 at the same redshift. 

For the AGN group at redshift 1.8-3.0 (z5-B), we could 
not find any evidence of a clustering signature. At this 
high redshift range, the number of galaxies which are 
brighter than the limiting magnitude is very poor, and the 
clustering signal easily disappears due to the foreground 
galaxies. 

The detection limit of galaxies for a lower redshift group 
is dimmer in absolute magnitude than that for a higher 
redshift group. From observations of the local Universe, 
it is known that brighter galaxies usually have stronger 
clustering. Zehavi et al. 2005 analyzed data from SDSS 
DR2 and derived the luminosity dependence of the galaxy 
clustering. They showed that the galaxy-galaxy correla- 
tion function increases continuously with luminosity, and 
above the characteristic luminosity (A/,. ^ —20.5) it 
increases more rapidly. The measured correlation lengths 
were 9.81± 0.39 h~^Mpc for Mr < —22, decreased rapidly 
down to 6.16±0.17 for —22 < Mr < —21, and decreased 
rather slowly down to 2.68±0.39 for -18 < Mr < -17. 
Thus it is expected that the clustering of galaxies around 
AGNs can be stronger at higher redshift due to the galaxy 
selection effect. 

In figure 20, the distributions of averaged number den- 
sity of bright (m < —21 + DM, left panel) and dim (m > 
—21 + DM, right panel) galaxies around the AGNs are 
shown for z4-B AGN group. The range of absolute mag- 
nitude for both samples can be inferred from the bot- 
tom panel of figure 18, which shows the distribution of 
m — DM{z) for R-band observations. They are expected 

to be in the range of —24 21 and —21 19 for the 

bright and dim galaxy samples, respectively. As can be 
seen from this figure, the excess number densities observed 
for dim and bright galaxy samples is almost comparable. 
The excess of source count at < 2.5Mpc is 9.3±1.2 
and 7. Oil. 2 for the dim and bright galaxy sample, re- 
spectively. Thus half of the clustering found for the z4-B 
sample is contributed from bright galaxies with absolute 
magnitude less than —21. 

It should be noted that this docs not necessarily mean 
that the clustering strength of dim galaxy sample and 



bright galaxy sample is equivalent, because the average 
numbers of detectable galaxies at the AGN redshift are 
different for the two samples. Since po for each of the 
bright and dim galaxy sample cannot be determined by 
the method as shown in section 3, we estimated the upper 
limit for the number density of bright galaxies (M < —21) 
by choosing smaller one of the followings: integral of the 
luminosity function up to M = —21 or po calculated for 
the whole galaxies as described in section 3. Then we 
obtained the upper limit of the average number density 
for the bright galaxies of z4-B sample as ^ 0.26 x 10~^ 
Mpc~^. The average density estimated for the whole 
galaxies sample is 1.3 x 10"'^ Mpc^^, thus the number 
density of the bright galaxies is less than 1 /4 of the dim 
galaxies (M > -21). 

The lower limit of correlation length for the bright 
galaxy sample was calculated as 18 h~^Mpc by adapting 
Po = 0.26 X lO"-^ Mpc~^ in the model fitting. This correla- 
tion length is significantly larger than that measured for a 
bright galaxy at the local Universe. Therefore the degree 
of increase in correlation length for z4-B relative to those 
for the low redshift samples is larger than that expected 
from the galaxy luminosity dependence observed in the 
local Universe. 

There would be another possibility that the redshift de- 
pendence is caused by selection bias to a type of galaxy, 
that is, early or late type galaxy. It is known that the 
early type galaxies are clustered stronger than the late 
type galaxies at the local Universe. Although we used 
UKIDSS IR data, the observations were too shallow to 
detect the high redshift galaxies and they were almost not 
used for the high redshift samples. Thus our high redshift 
dataset is biased to blue-colored late type galaxies. This 
is contrary to the expectation that the larger correlation 
length for z4-B group is due to the selection bias to the 
redder galaxies which are strongly clustered in the local 
Universe. 

6. Discussion 

We showed that the low redshift AGNs (z = 0.3-0.9, zl 
and z2) had clustering properties similar to that of a typ- 
ical local galaxy. This would imply that these AGNs are 
located in dark matter halos which have a mass similar 
to that occupied by typical galaxies at the local Universe, 
provided that there is a correlation between the cluster- 
ing property and the mass of the dark matter halo as 
suggested by many numerical simulations. 

No significant AGN luminosity dependence was ob- 
served for z2 and z3 redshift samples. Since the luminosity 
of AGN is not constant and evolving with time, it may be 
possible for AGNs with different luminosities to have sim- 
ilar environments and to be powered by the same mecha- 
nism at some epoch of their evolution. The QSO lifetime, 
tq, which is a timescale for the most luminous phase of 
the activity, is expected to be in the range from 10^ - 10^ 
yr (Martini 2004). For a duration much longer than the 
period of the most luminous phase, AGNs arc obscured 
by surrounding gas and dust, and thus they are observed 
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as low luminosity AGNs or other types of galaxies such as 
sub-mm galaxies, or they cannot be observed completely. 
Therefore, no detection of the luminosity dependence of 
the clustering supports the theoretical expectation of the 
short QSO lifetime. 

We also showed that the bright AGNs at redshift 1.3-1.8 
(z4-B) have relatively denser environments than typical 
local galaxies, and that the clustering strength was much 
larger for bright galaxies {M < —21) than dim galaxies 
(M > — 21). The lower limit of the cross correlation length 
of AGN and the bright galaxies was 18 h~^Mpc, which 
is extremely large compared to the local bright galaxies 
measured by Zehavi ct al. (2005). Considering that our 
dataset is biased to blue star-forming galaxies at higher 
redshift, we deduce that the clustering of blue star-forming 
galaxies is enhanced around the high redshift AGNs. 

Under these circumstances, the chance of a major 
merger is expected to be enhanced relative to that at low 
redshift, thus the major merger could be a preferred mech- 
anism to trigger the AGN activity at this redshift range. 
This result indicate that there exists redshift dependence 
for the environment of AGNs, that is, at higher redshift 
(z >~1.3) larger fraction of AGNs are produced in a more 
crowded region with blue star-forming galaxies, while at 
lower redshift they are mostly produced in an environment 
which is common at the local Universe. AGNs at low red- 
shift might be powered by different mechanism from high 
redshift AGNs, e.g. a secular process such as bars or other 
disk instabilities. 

The w(rp) profiles of z4-B show deviations from a power 
law function, which shows a flat over density at Vp <2.5 
Mpc. This property could be attributed to the distribu- 
tion of AGN location in a dark matter halo; there may be 
a considerable number of AGNs which have an offset of 
^2.5 Mpc from the center of galaxy distribution or a dark 
matter halo. 

The w(?'p) for bright AGN samples have a tendency to 
show a deficit of galaxy at the smallest distance bin. This 
feature is not seen for any of the dim AGN samples. This 
also might be explained by distribution of AGN location 
inside the dark matter halo. Galaxy formation near the 
bright AGN may also be suppressed due to the strong 
radiation from the AGN. 

In this work, the error of ro is large compared to the ex- 
periments measuring the redshift of all the target galaxies. 
This is mostly due to the uncertainty of determination po, 
thus the accuracy of po needs to be improve. For that pur- 
pose, we need to calibrate the po estimated as described 
in section 3 experimentally. Once the po estimation is cal- 
ibrated and the accuracy is evaluated for each AGN field, 
it would be possible to reduce the uncertainty of po by a 
factor two. 

7. Conclusions 

Using the Japanese Virtual Observatory, we were able 
to measure the clustering property for ~ 1.809 AGNs at 
redshifts 0.3-3.0 with optical brightness of My = —30 to 
—20. The AGN samples were divided into five redshift 



groups (zl-z5); and each redshift group was further di- 
vided into two brightness groups (e.g. zl-D, z2-D, z2-B 
...). 

We found evidence of clustering except for the highest 
redshift group (z5-B). The correlation lengths between the 
AGNs and galaxies were measured for each group. We 
found that the clustering strength of the high redshift 
samples (z4-B) was larger than that of the low redshift 
sample (zl-D), and it was nearly consistent with that of 
QSO clustering at the same redshift. We also showed that 
bright galaxies of absolute magnitude {Mab < ~21) were 
more densely clustered than the dimmer galaxies for the 
case of z4-B AGN group. Considering that, at higher red- 
shift (z > 1), the data used in this work are based on 
optical observations and blue galaxies are expected to be 
preferentially detected in our dataset, we can deduce that 
blue star-forming galaxies are highly clustered around the 
AGNs at these redshifts. This clustering feature is differ- 
ent from that seen for the low redshift AGNs. 

The deviation of a; (rp) from the power law function was 
seen in the result of z4-B samples, and there was an indi- 
cation for suppression of cross-correlation at the smallest 
distance bin for the bright AGN samples. These features 
could be related to the distribution of AGN location in a 
dark matter halo, and could be related to fueling mech- 
anism of the AGN. If a significant fraction of AGNs was 
located offset from the center of a dark matter halo or a 
galaxy distribution, the deviation from power law could 
be seen. It is also possible that the galaxy formation near 
the AGN was suppressed due to the strong radiation of 
the AGN. 

We showed that the analysis method presented in this 
paper can be a powerful probe for the small scale environ- 
ment of the high redshift AGNs. Our method can be used 
for tracing dim and blue star-forming galaxies at small 
scale (r < 1 Mpc), while the spectroscopic observations 
provide a precise information about large scale clustering 
properties for bright and red early type galaxies. These 
two types of observations, therefore, are complementary, 
and can provide important clues to understanding the ori- 
gin of the AGN activity. 

To derive more definite evidence for the AGN-galaxy 
clustering with this method, more observational data are 
needed in view of the quantity and quality. Highly homo- 
geneous observations of a huge number of AGN fields will 
be conducted with the Hyper Suprime-Cam, which will be 
installed on the Subaru Telescope around 2011. This sur- 
vey project will give more definite observational evidence 
on the evolution of an AGN. 
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Table 1. Distance ranges for background estimation in the calculation of Bqq 



redshift group 


redshift 


range 


^bg,min 


^bg,max 








Mpc 


Mpc 


zl 


0.3 - 


0.6 


1.5 


2.5 


z2 


0.6 - 


0.9 


2.0 


3.5 


z3 


0.9 - 


1.3 


3.0 


5.0 


z4 


1.3 - 


1.8 


4.0 


6.0 


z5 


1.8 - 


3.0 


4.0 


6.0 



Table 2. Number of AGNs selected from SDSS, 2dF, other optical, X-ray and radio observations for each redshift and AGN 
luminosity group. 



group 


SDSS * 


2dF t 


UV-OPT-IR ^ 


XRAY 


RADIO II 


Total 


zl-D 


478 (73.5%) 


51 (7.8%) 


96 (14.7%) 


20 (3.1%) 


6 (0.9%) 


651 


z2-D 


334 (79.5%) 


66 (15.7%) 


13 (3.1%) 


7 (1.7%) 


(0.0%) 


420 


z3-D 


68 (68.0%) 


20 (20.0%) 


9 (9.0%) 


3 (3.0%) 


(0.0%) 


100 


z2-B 


101 (89.4%) 


4 (3.5%) 


2 (1.8%) 


2 (1.8%) 


4 (3.5%) 


113 


z3-B 


234 (85.4%) 


32 (11.7%) 


4 (1.5%) 


2 (0.7%) 


2 (0.7%) 


274 


z4-B 


107 (75.4%) 


21 (14.8%) 


10 (7.0%) 


4 (2.8%) 


(0.0%) 


142 


z5-B 


76 (69.7%) 


21 (19.3%) 


8 (7.3%) 


4 (3.7%) 


(0.0%) 


109 


total 


1398 (77.3%) 


215 (11.9%) 


142 (7.8%) 


42 (2.3%) 


12 (0.7%) 


1809 



• Number of AGN samples contained in the SDSS DR5 QSO catalog (4th Ed.) (Schneider et al. 2007). This number 
includes samples contained also in the 2dF QSO catalog. 

t Number of AGN samples contained in the 2dF QSO catalog (Groom ct al. 2004) and not contained in the SDSS QSO 
catalog. 

f Number of AGN samples selected based on the optical properties. References of each AGN are: Zakamska et al. 2003, 
Reyes et al. 2008, Abazajian et al. 2004, GoUinge et al. 2005, Schneider et al. 2005, Prcscott ct al. 2006, Barton et al. 2006, 
Zhan and Chen 1989, Williams et al. 2002, Richards et al. 2001, Hall et al. 2000, Veron et al. 1990, Rowan-Robinson et al. 
1990, Mcintosh et al. 2004, Liu et al. 1999, Leipski et al. 2005, La Franca et al. 1992, Kniazev et al. 2004, Hewett et al. 
1995, Cristiani et al. 1995, Boyle et al. 1990, Hennawi et al. 2006, Hall et al. 1996, Zhan et al. 1987, Schneider et al. 1994, 
Hewett et al. 1991, Sharp et al. 2002, Abazajian et al. 2009, Hoag et al. 1982, Usher et al. 1983, Monk et al. 1988, He et al. 
1984, Dressier & Gunn 1992, Cristiani et al. 1989 Zhou et al. 2006 Schmidt et al. 1986 

§ Number of AGN samples selected based on X-ray observations. References of each AGN are: Silverman et al. 2005, 
Steffen et al. 2004, Lehmann et al. 2000, Watanabe et al. 2002, Schmidt et al. 1998, Molthagen et al. 1997, Lehmann et al. 
2001, La Franca et al. 2002, Hasinger et al. 2002, Grindlay 1993, Gioia et al. 2003, Eckart et al. 2006, Dewangan et al. 
2001, Crawford et al. 2002, Caccianiga et al. 2004, Branduardi-Raymont et al. 1985, Bauer et al. 2000, Barger et al. 2002, 
McHardy et al. 1998, Miller et al. 2002, Georgantopoulos et al. 2004, Fiore et al. 2003, Bower et al. 1996, Bade et al. 1995, 
Mason et al. 2000, Arp et al. 2001, Grampton et al. 1997, 

II Number of AGN samples selected from Radio observations. References of each AGN are: Becker et al. 2001, Glikman 
et al. 2004, Healey et al. 2008, Hes et al. 1996, Lynds et al. 1966, Marziani et al. 1996, Rowan-Robinson et al. 2004, White 
et al. 2000, Wold et al. 2000 Buchalter et al. 1998 



Table 3. Number of AGN samples for each AGN group and observation band. 



group 


B 


V 


R 


I 


i' 


z' 


K 


zl-D 


13 


7 


21 


9 


13 


4 


584 


z2-D 





6 


18 


16 


4 


2 


374 


z3-D 


2 


3 


4 


6 


5 


4 


76 


z2-B 


1 


1 


9 


7 


3 


1 


91 


z3-B 


10 


15 


25 


19 


11 


8 


186 


z4-B 


11 


15 


45 


32 


21 


9 


9 


z5-B 


12 


19 


38 


18 


12 


5 


5 
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Table 4. Fitting parameters of projected cross correlation analysis 



group 




z 


z* 


My 


t 


"AGN* 

mag 


I'o" 
h^^Mpc 


Mpc"2 


Pq* 
Mpc^^ 


zl-D 


0.3 


- 0.6 


0.46 


-25.0 - 


-20.0 


651 


^- ' -0.7 


22.1 ± 0.07 


4.3+'^^ 


z2-D 


0.6 


- 0.9 


0.73 


-25.0 - 


-20.0 


420 


5.8ti? 


10.1 ± 0.05 




z3-D 


0.9 


- 1.3 


1.01 


-25.0 - 


-20.0 


100 




11.2 ± 0.07 


^■"-0.5 


z2-B 


0.6 


- 0.9 


0.79 


-30.0 - 


-25.0 


113 




13.9 ± 0.10 


9 q+0.4 
^•"^-0.6 


z3-B 


0.9 


- 1.3 


1.06 


-30.0 - 


-25.0 


274 


5 1+2-7 
3-1^-1.5 


11.1 ± 0.04 


1 4+0-4 


z4-B 


1.3 


- 1.8 


1.54 


-30.0 - 


-25.0 


142 


11 1+^-1 
ll.i_2 7 


19.4 ± 0.07 


1 q+0.4 
^■^-0.6 


z5-B 


1.8 


- 3.0 


2.07 


-30.0 - 


-25.0 


109 


< 13.1 


14.4 ± 0.07 


4Q+0-19 
"■+-'-0.24 


zl-D (IR) 


0.3 


- 0.6 


0.46 


-25.0 - 


-20.0 


584 


6.8±l:|i 


7.73±0.05 


9 9+0.4 
^•^-0.6 


zl-D (OPT) 


0.3 


- 0.6 


0.46 


-25.0 - 


-20.0 


67 


3.0+1-^ 


147.0±0.67 


22.3^2-^ 


z2-D (IR) 


0.6 


- 0.9 


0.73 


-25.0 - 


-20.0 


374 




3.50±0.03 


0.56l°i? 


z2-D (OPT) 


0.6 


- 0.9 


0.74 


-25.0 - 


-20.0 


46 


5 0+1-^ 


63.5±0.34 





t K-corrected absolute magnitude range 
t number of AGNs 

il correlation length, the quoted error contains systematic error due to uncertainty of po and one sigma statistical error. 

§ projected galaxy number density for background 

* average galaxy number density at the AGN redshift 
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Fig. 1. The distribution of the average number density po of detectable galaxies at the AGN redshift z. Data of po > lO"'^ Mpc~^ 
among the 2689 samples are plotted. The open circles and crosses represent AGN samples for which the galaxy sample is derived 
from the UKIDSS and the Suprime-Cam data, respectively. The horizontal dashed line represents the threshold for data selection. 
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Fig. 2. The logarithmic distribution of clustering coefficient Bqg for 2023 samples which are selected with the criterion po > 10~^ 
Mpc~'^. The top panel is for negative Bqq and the bottom panel is for positive Bqg- The vertical dashed lines represent the upper 
limits for data selection. 
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Fig. 3. The distribution of reduced of number density at radial distance from 1 Mpc to rj^g j^^jj^ for 1976 AGN samples. The 
vertical dashed line represent the maximum value for sample selection. 
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Fig. 4. The distribution of maximum deviation of number density measured in the unit of a standard deviation at radial distance 
from 1 Mpc to r-^^^ j^^jj^ for 1976 AGN samples. The vertical dashed line represent the maximum value for sample selection. 
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Fig. 5. K-corrected V band absolute magnitude vs redshift of the AGNs used in this work. Open circles represent AGN samples 
for which the galaxy sample is derived from the UKIDSS data, and the crosses represent AGN samples for which the galaxy sample 
is derived from the Suprime-Cam data. 
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Fig. 6. Redshift distributions for dim AGN samples (top) and bright AGN samples (bottom). The dashed histogram represents 
AGN samples for which the galaxy sample is derived from the UKIDSS data, and the solid line histogram represents AGN samples 
for which the galaxy sample is derived from the Suprime-Cam data. 
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Fig. 7. The distribution of 1,809 AGNs analyzed in this work. The coordinates are in equatorial. 
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Fig. 8. The distribution of AGN coordinates. The top panel is for the right ascension of AGN, the left bottom panel is for the 
declination, and the right bottom panel is a close up of the declination distribution, showing detailed distribution of low counts. 
Open histogram is for all the AGNs and the shaded histogram is for AGNs for which the Suprime-Cam data was used. 
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Fig. 9. Examples of the observed magnitude distributions and their model functions. The left histogram shows the UKIDSS K-band 
data and the right one shows the Suprime-Cam R-band data. 
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Fig. 10. Distribution of nijijjjij and AGN redshift. The open circles correspond to AGNs for which the galaxy sample is derived 
from the UKIDSS data, and the crosses correspond to AGNs for which the galaxy sample is derived from the Suprime-Cam data. 
The dashed line represents apparent magnitude which corresponds to absolute magnitude of —21 at the redshift. 
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Fig. 11. The distributions of number density of detected sources and tlie projected correlation functions for dim AGN groups. The 
left panels shows the averaged galaxy number density as a function of projected distance from the AGN, and the corresponding 
projected correlation functions are shown in the right panels. The dashed lines in the left panels represent fitting parameter njjg. 
The solid lines represent the model function fitted to the observational data. 
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Fig. 12. The distributions of number density of detected sources and the projected correlation functions for bright AGN groups. 
The other explanation of this figure is the same as figure 11. 
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Fig. 13. The logarithmic distribution of clustering coefficient Bqq for 142 z4-B samples. The top panel is for negative Bqg and 
the bottom panel is for positive Bqq. 



28 



Y.Shirasaki et al 



[Vol. , 



24 



23 - 



a. 21 



20 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

24-B 



(IBq^I <= 1000) J 




^C) 1 1 I I I I I I I I 1 1 I I I I I I I I 1 1 I I I I I I I I 1 1 I I I I I I I I 1 1 I I I I I I I I 1 1 I I I I I I I I 

1 2 3 4 5 6 

fp [Mpc] 



Fig. 14. Averaged number density around z4-B AGNs which satisfy the condition |BqgI ^ 1000. 
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Fig. 15. Averaged number density around z4-B AGNs not located in the Sextans field (left) and AGNs located in the Sextans field 
(right). 
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Fig. 16. Measured correlation lengtli, rp, as a function of redshift. The error bars of our results indicate sum of the systematic 
error estimated from uncertainty of po and the statistical error of one sigma. The closed triangles represent the result for dim AGN 
groups, and the closed circles represent the result for bright AGN groups. The result for z5-B are shown with the upper limits. 
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Fig. 17. The its of the galaxies around AGNs (cg aGN-G)i QSOs (eg q), and galaxies (cg (-;) are shown as a function of redshift. 
The error bars of our results indicate sum of the systematic error estimated from uncertainty of po and the statistical error of one 
sigma. 



No. ] AGN and Galaxy Clustering at z = 0.3 to 3.0 31 



c 

o 
O 



10^ - 



n I I I I I I I I I I I I I I I I I I I I I I r 
z1-D 

Suprime-Cam R band 



10^ - 



10^ ^ 



10' - 



- UKIDSS K band 




10^ 



r 1 l| 



iiipi 

I II if iiii 

nil .1 I I I I I L 



I 



I I' 

I'li'i',," 

I I I I i'!"l'!'i.'i 



-40 -35 -30 -25 

m - DM(z) 



-20 



-15 



c 

o 
O 




-30 -25 
m - DM(z) 



Fig. 18. The distribution of m — DM{z), where m is an apparent magnitude of detected objects, and DM{z) is a distance modulus 
corresponding to the AGN rcdshift. The top panel shows the distribution for zl-D galaxy samples observed with Suprime-Cam in R 
band (solid line) and in UKIDSS K band (dashed line). The bottom panel shows the distribution for z4-B galaxy samples observed 
with Suprime-Cam in R band. 
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Fig. 19. Averaged number density around z3-B and z4-B AGNs. 
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Fig. 20. Averaged number density of dim sources [m > —21 -f DAI) (left panels) and of bright sources (m < —21 H- DM) (right 
panels) around z4-B AGNs, where DM represents a distance modulus corresponding to AGN redshift. The solid lines represent 
model function of equation (13) fitted to the observation data. The fit was done for distance range of 0.1-2.0 Mpc for the case of 
dim sources, while 0.0—2.0 Mpc for bright sources. The horizontal lines represent the fitting parameter n^^g. 



